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The validity of the affine assumption in model flexible polymer networks is explored. To 

T-H 

^ this end, the displacements of fluorescent tracer beads embedded in polyacrylamide gels are 

t-H 

O quantified by confocal microscopy under shear deformation, and the deviations of these dis- 

T~H 

^ placements from affine responses are recorded. Non-affinity within the gels is quantified as a 

t-H 

t-h function of polymer chain density and cross-link concentration. Observations are in qualita- 

o 

tive agreement with current theories of polymer network non-affinity. The measured degree 
• ^ of non-affinity in the polyacrylamide gels suggests the presence of structural inhomogeneities 

X 

which likely result from heterogeneous reaction kinetics during gel preparation. In addition, 
the macroscopic elasticity of the polyacrylamide gels is confirmed to behave in accordance 
with standard models of flexible polymer network elasticity. 
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1 Introduction 



Affine deformation is an essential assumption in the classical theory of elasticity. In the classical 
theory, deformation is assumed to be distributed homogeneously in the material so that strain is 
spatially constant at all length scales. The affine assumption permits elastic properties of cross- 
linked polymer networks to be readily derived from theories of rubber elasticity based on the 
entropy of a single polymer chain in the network. In practice, however, such affine deformations 
only occur in perfect crystals under very small deformation. In polymer networks, especially net- 
works composed of semi-flexible or rigid filaments, the microscopic network deformations should 
be non-affine below certain length scale. 

Non- affinity can arise from different sources. In near- ideal flexible polymer melts, deforma- 
tions might be expected to be affine on length scales much larger than the average mesh size and 
non-affine at lengths scales of the order of the mesh size or smaller. 1 Random thermal fluctuations 
of the cross-link junctions, along with thermal undulations of the polymer chains may also lead to 
non-affine behavior in polymer gels. Inhomogeneities introduced into the network micro-structure 
during sample preparation can also introduce non-affine responses; such inhomogeneities might 
be expected to be a function of reaction kinetics and other sample preparation parameters. 2 

Over the years, the connection between shear deformation and non- affinity has been explored 
theoretically in a wide range of materials including rubber-like spatially homogeneous elastic me- 
dia,^ entangled or cross-linked polymer networks,'^' semi-flexible polymer networks with rigid^ 
and flexible cross-links,^ stiff rod networks,^ biopolymers,^ amorphous systems 12 and foams.^ 
Indeed, it has been proposed that non-linear elasticity in polymer networks has its origin in non- 
affine responses.^ In spite of continued interest in this problem and its fundamental importance, 
relatively little experimental quantification of the non-affine phenomenon has been carried out in 
semi-flexible biopolymer gelsj 15 * 20 ^ and we are not aware of any non-affinity studies for the simple 
case of flexible polymers. Experiments along these lines should provide benchmarks for future 
understanding of the subject. 

This paper describes an investigation of non-affine shear deformations in a model flexible poly- 
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mer gel: polyacrylamide gels with bisacrylamide cross-links. Polyacrylamide is well suited for the 
investigation because it is comparatively well-controlled, and its stiffness is tunable by the number 
of bisacrylamide cross-links. As part of this study, macroscopic rheological measurements are car- 
ried out to confirm the simple rubber-like elastic character of these networks. Then deformation 
fields in the gels under external shear stress are characterized by measuring the displacements of 
fluorescent beads entrapped in the gels. Bulk rheology and confocal microscopy are used in con- 
cert for the latter study. A non-affine parameter, A, is defined to quantify the degree of non-affinity 
in the displacement field. A is measured as a function of bead size, polymer chain density, and 
cross-link density in the gels. We test simple predictions of a recently developed theory of non- 
affinity in random elastic media 2 and obtain estimates for the fluctuations in elastic modulus of the 
gels from A. 

2 Experiment 

2.1 Sample Preparation 

The polyacrylamide (PA) gel is prepared by polymerizing acrylamide monomers and bisacry- 
lamide (bis) cross-links in aqueous 50 mM HEPES buffer at pH = 8.2, using free-radical polymer- 
ization reaction initiated by 0.1% weight/weight (w/w) ammonium persulphate (APS) and 0.3% 
w/w N, N, N', N'- tetramethylethylenediamine (TEMED). (Here the percent of X w/w equals the 
mass in grams of X per 100 grams of solution.) Fluorescent polystyrene tracer beads are mixed 
into our solution at a concentration of 0.004% weight per volume (w/v), before the addition of 
bis cross-links. (Here the percent of X w/v equals the mass in grams of X dissolved/suspended 
in 100 milliliters of solvent.) Thus, a tracer bead concentration of 0.004% w/v is attained by dis- 
solving 0.004 gram of tracer beads in 100 milliliters of water. This procedure helps to distribute 
the beads uniformly throughout the polymer network. Internally labeled and carboxylate-modified 
fluorescent polystyrene micro-spheres of various diameters are used for this purpose, viz., 0.6 /im, 
1 /im (Molecular Probes, California, USA), and 1.5 /im (Bangs laboratories Inc., Indiana, USA). 

3 



Acrylamide (7.5%, 15% w/v and bisacrylamide (0.03 - 0.12%w/v) concentrations are systemati- 
cally changed to study the effects of polymer concentration, cross-link density and mesh size on 
the polymer network rheology. 



2.2 Rheology 

Rheology measurements are performed using a stress-controlled Bohlin Gemini rheometer (Malvern 
Instruments, UK), with a cone and plate geometry of 4° cone angle and 20 mm diameter. Sam- 
ples are prepared in situ so that good contact is routinely established between the sample sur- 
faces and the rheometer plates to prevent slippage at high strains. The shear modulus (G ; ) and 
loss modulus (G") for each sample during the process of polymerization are monitored using low 
strain amplitude (yo= 0.01) and low frequency (/ = 0.1 Hz) oscillatory shear measurements. The 
polymerization reaction proceeds for ~ 30 minutes, with the elastic and viscous moduli attaining 
steady-state values in less than 10 minutes. Care is taken to prevent solvent evaporation by sealing 
off the sample from the sides with a low density, low viscosity (~ 50 mPa-s) silicone oil. The 
elastic and viscous moduli, G' and G", respectively, for these gels are measured as functions of 
frequency, amplitude and temperature. These measurements are intended to confirm that the gels 
behave in accordance with the existing theories of flexible polymer networks. 11 ^ 1 A set of control 
experiments are performed on the PA gels, with and without the tracer beads, to further confirm 
that macroscopic properties of the gels are not altered by the addition of the tracer beads. 



2.3 Confocal Microscopy 

Microscopic deformation of the PA gels under shear is studied by tracking tracer bead displace- 
ments in the sample using confocal microscopy. A VTeye confocal system (VisiTech International, 
UK) is used in conjunction with an inverted Eclipse TE200 microscope (Nikon Instruments, USA) 
for this purpose. The lower plate of the rheometer is replaced by a home -built transparent sam- 
ple holder and is mounted on the microscope to permit visualization of the samples under shear 



[Fig. [figure][l][]i;a)] 




Figure 1: Experimental Setup, (a) Experimental schematic, (b) Sketch of the non-affine displace- 
ments of tracer beads. (xoi,yoi,Zoi) and (xi,yi,Zi) mark the positions of a tracer bead without and 
under shear, respectively. Dashed arrow indicates affine displacement, d a i, of tracer bead in the 
direction of shear (x-axis). di is the measured displacement of the tracer bead. u X i, u y u u z t indi- 
cate the non-affine deviations along the x, y and z axes respectively, the non-affine 
deviation. 
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A 60 x water objective (NA = 1.2) is used to visualize the sample over a depth of 100 /im. 
3D stacks map the entrapped tracer beads (70 /im x 70 /im x 60 /im) using the confocal setup 
with and without applied shear. The step size of the 3D stacks is varied from 100 nm to 200 nm 
for tracer bead sizes ranging from 600 nm to 1.5 /im. A wide range of shear strain, up to 50% 
amplitude, is applied. 

The image stacks are processed using relatively standard Matlab routines, which determine 
the beads' positions with subpixel accuracy,^.^ For each stress value, a set of two image stacks 
are taken, one with shear and one without. The 3D locations of beads in a gel without external 
shear stress are determined as (xoi,yoi,Zoi), for i = l,2,...,N, where N is the number of tracked 
beads. The centroids of the corresponding /V beads in the image stack under shear stress are 
measured too, as ix^y^Zi); for convenience the direction of shear is taken to be along the x axis 
in the figure. The displacements of tracer beads are then calculated from the tracking results as 
di = (xi—xouyi —you Zi — zoi). The system permits the displacements of tracer beads to be measured 
with a spatial resolution of 50 nm. On average, 30 beads are tracked in each 3D stack. 



2.4 Non-affine Parameter 

A measure of the degree of non-affinity is provided by the non-affine parameter A, which is defined 
in Ref.^as: 

1 N 

Here Uj = di — d a i, is the deviation of the measured tracer-bead displacement, di, from the affine 



displacement, d a i [Fig. | [figure] [ 1 ] [] 1 [ b)] . 

For a perfect shear deformation along the x-axis, the affine displacement d a i would be in the 
direction of shear only- the y and z components must be zero. We measure the resultant strains 
along all three component axes, y x , y y and y z by fitting the x, y and z components of di to linear 
functions of zqu as seen from a sample PA gel (7.5% acrylamide, and 0.03% bis) under an applied 
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Figure 2: (a) Experimentally measured displacements of tracer beads in the direction of shear, 
di, that has been decomposed along x-(red circles), y-(blue crosses) and z- (green pluses) axes, as 
a function of the distance, zqi from the fixed lower plate of the rheometer. This sample is 7.5% 
acrylamide and 0.03% bis PA gel. The solid lines give the strains, y x , y y , and y z obtained from 
their linear fits. Note that y y , and y z are w 0. (b) The distribution of non-affine deviations of tracer 
beads for the same sample PA gel shown in Fig. |[figure][2][]2f a), at y = 0.3, decomposed along 
the x-, y- and z- axes. The measurements are normally distributed around the affine displacement 
position, as indicated by the solid curves. 
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strain of 7 = 0.3 in Fig. [figure] [2] [] 2 a). The real strain on the sample is determined as 7 = 
Jrl + Yy+Yz- The x, y and z components of the affine displacement vector, d a i are then calculated 
as zoijx, ZQiJy, and zwYz- Note that the y and z components, both perpendicular to the direction of 



shear, do not vary as a function of zu resulting in y y and y z w . Fig. [figure] [2] [] 2 b) plots the 
distribution of non-affine deviations, u x , u y , and u z for the same sample gel, along the x, y, z axes 
respectively, for the same strain of 7 = 0.3 as seen in Fig. [figure] [2] [] 2 a). |w ; -| 2 is calculated as 
(xi - x 0i - jxzoi) 2 + (yi - yoi - 7-zo/) 2 + (zi - zoi - y z zoi) 2 . 

The non-affine parameter, A is then defined in terms of these variables as 

1 N 

A = m £ K x »" - X 0i ~ Y* Z 0i) 2 + ~ ~ ^0/) 2 + (Zi ~ Z 0i ~ YzZOif] • (1) 

3 Results 

3.1 Bulk Rheology Measurements 

The PA gels used in our experiments are solid-like materials with G' ranging from 7.7 x 10 2 Pa 



to 1.5 x 10 4 Pa, 2 to 3 orders of magnitude larger than the G" . In Fig. [figure] [3] []3 a), G' and 
G" of a gel made of 7.5% acrylamide and 0.06% bisacrylamide are plotted as functions of the 
amplitude of the oscillatory shear strain at oscillation frequency, / = 0. 1 Hz. G' is approximately 
100 times larger than G" . Moreover, both G' and G" are independent of the applied shear strain 
for strains up 7 = 0.5, confirming the linear elastic response of PA gels. The frequency response 
of PA gels is characterized by measuring G' and G" at oscillatory strains with amplitude 70 = 0.01 
and frequency ranging from 0.1 Hz to 100 Hz. Within this frequency range, G' remains constant, 
and G" increases with increasing frequency (data not shown). 

The elastic moduli of our PA gels vary linearly with bisacrylamide concentration and sample 



temperature. Cross-link and monomer concentration trends are shown in Fig. [figure] [3] [] 3 1>). 
Notice, when the bisacrylamide concentration increases from 0.03% to 0.12%, G' for gels with 
7.5% acrylamide increases linearly from 7.7 x 10 2 Pa to 4.9 x 10 3 Pa. Similarly, G' for 15% 
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Figure 3: Rheology of polyacrylamide gels, (a) G' of sample PA gel is two orders of magnitude 
larger than the G" and theses values remain constant over a wide range of applied strain. Data 
are shown for a gel with 7.5% acrylamide at 0.03% bisacrylamide cross-link concentration, at an 
oscillatory frequency of 0.1 Hz. (b) G' of 7.5% and 15% polyacrylamide gels as a function of 
cross-link concentrations. Error bars denote standard deviations which are less than 2% of the 
mean elastic moduli. The solid lines indicate linear fits to the data. Note that the overall moduli of 
the gels with 7.5% acrylamide are significantly lower than that of 15% acrylamide for comparable 
cross-link density, (c) G' as a function of temperature (red line is the linear fit). Data are shown for 
a PA gel with 7.5% acrylamide with 0.09% bisacrylamide. 

acrylamide PA gels increases from 1.6 x 10 3 Pa at 0.005% bis to 1.5 x 10 4 Pa at 0.05% bis con- 
centration. We also investigated the temperature dependence of the network elasticity within the 



attainable temperature range of the rheometer, i.e. 5°C < T < 90°C. In Fig. [figure] [3] []3 c), we 
show that G' from the gel made of 7.5% acrylamide and 0.09% bisacrylamide increases linearly 
with sample temperature. This linear dependence of G' on cross-link concentration and sample 
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temperature follows the predictions of classical rubber elasticity theory. Note that the slopes of the 
linear fits of G' as a function of cross-link concentrations for the 7.5% acrylamide is lower than 
that of 15% acrylamide PA gels. We suggest that some bis molecules form efficient cross-links 
and others do not, and that this difference in the slopes of G' versus cross-link concentrations for 
7.5% and 15% acrylamide is due to the difference in effectiveness of the bis molecules in forming 
efficient cross-links, which increases with increasing monomer concentration. We discuss these 



effects further in Section [subsection] [ 1 ] [4]4. 1 



3.2 Non-affine parameter, A scales as the square of the applied strain 

Confocal microscopy is used to visualize and record the displacements of the fluorescent tracer 
beads entrapped within a (70 fim x 70 /im x 60 /im) volume in the PA gel. Since the tracer beads' 
size of ~ 1 \lm is much larger than the average mesh size of the PA gel, free Brownian motion is 
suppressed. Within this small volume, located approximately 1 cm from the axis of rotation, the 
macroscopic shear strain applied to the beads can be approximated to be unidirectional. 

In Fig. |[figure][2][]2f a), bead displacements along the x, y, and z axis, are plotted as a func- 



tion of zoi, the distance between the beads and the bottom surface. The displacements along the 
direction of shear, viz., the x-axis, increase linearly with zot, as expected from macroscopic shear 
deformation. Fitting d x to a linear function of zoi yields the strain y x « y. d y and d z , both perpen- 



dicular to the shear direction, are independent of zoi as shown in Fig. [figure] [2] [] 2 'a). Also notice 



from Fig. [figure] [2] [] 2 ;b) that the non-affine displacements along each axis, viz., u x , u y and u z , 
are much larger than the resolution of our system in the xy— plane (~ 50 nm) and comparable to 
that along the z axis (~ 80 nm) and are normally distributed with mean value zero, i.e., distributed 
around the affine displacement positions. These uncertainties in tracking lead to a noise floor in 
A ~ 0.007 \im 2 . 

The non-affine parameter A is readily computed from the measured bead displacements using 
Eq. ([T]) for PA gels (7.5% and 15% acrylamide and a range of bis acrylamide concentrations). In 



Fig. [figure][4][]4^a), A increases with applied strain y, and clearly scales as y . 
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Figure 4: (a) The non-affine parameter scales as the square of the external strain field, as seen 
for sample polyacrylamide gels at 7.5% acrylamide and 0.03% bis (top), and 15% acrylamide and 
0.01% bis (bottom). The dashed lines give best linear fits to the data. Error bars represent standard 
deviation of measurements of strain and non-affinity, the latter being smaller than symbol size, (b) 
Strain-normalized non-affine parameter, 4 for sample PA gels at 7.5% and 15% acrylamide are 
plotted at varying bis concentrations. The data points and error bars represent the average value 
and standard error of measurements from different samples prepared in the same manner. The 
dashed line indicates the average 4 calculated from all data points in the figure. 

In Ref., DiDonna and Lubensky developed a perturbation theory for non-affine deformations 
in solids with random, spatially inhomogeneous elastic moduli. They characterized the non-affine 
deformations using the non-affinity correlation function 



<3ij{x,x') = (Ui(x)uj(x')), 



(2) 



where u(x) is the nonaffine displacement field, and (■} represents the average over randomness 
in the elastic moduli (i.e., a disorder average). Because the disorder averaged quantities are 
translationally and rotationally invariant in the gel we consider, the correlation function only 
depends on the distance \x — and thus it is characterized by the Fourier transform Q{q) 
j d(x — x')Gu(x,x')e~ iq '( x ~ x 'K In Ref. 2 it is proved that this correlation function is related to the 
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correlation function of the inhomogeneous elastic modulus K as 

Giq) ~ (3) 

where K is the disorder averaged elastic modulus, and A K (q) is the Fourier transform of the 
spatial correlation function of the elastic modulus K. 

In this theory, the zeroth order problem concerns elastic deformations in a homogeneous media 
of elastic modulus K, and the randomness in K is treated as a perturbation from this homogeneous 
state. To first order, the driving forces of the non-affine deformations are thus proportional to the 
zeroth order deformations, which are proportional to y. Therefore, to first order in perturbation 
theory, Q{q) is proportional to y 2 . 

The non-affine parameter A defined in the present experiment corresponds to Qu(x, x), 



a pi \ f d 3 q n( f d 3 q fA K (q) 



(4) 



It is clear from this equation that A y 1 . In our experiment the relation A y 2 is verified, as shown 



in Fig. [figure][4][]4^a). This fairly robust relation has also been found in non-affine correlation 
functions of, for example, flexible polymer networks^ and semiflexible polymer networks at small 
strain.™ 

The quantity 4, which is independent of strain y, provides a good measure of the degree of 



non-affinity of the sample. We shall refer to this quantity, 4* as me strain-normalized non-affine 
parameter. 4 is calculated for each sample as follows: For a particular strain, A is calculated by av- 



eraging the square of the non-affine displacements, uf, for all tracer beads in the sample; typically 
we carried out multiple shear measurements at the same strain (see Section [subsection] [3] [4]4~3] ), 



and the displacement data from all particles in all repeated shear measurements were averaged 
together to derive the mean A and its standard deviation. The resultant 4 data were then fit to a 



A 
f- 

comparable to the noise floor of the measurements in A. Standard deviations for the slopes were 



linear function. The slope of the linear fit gives 4 for the sample; the intercept form the fitting is 
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also derived. We use this parameter, A versus y, which represents an intrinsic material property, 
for comparisons among samples prepared at different times or under different conditions. 

The 4 value calculated for each sample, along with its constituent acrylamide and bis concen- 
tration is listed in Table 2 in the supporting material available online. Fig. |[figure] [4] []4{ b) plots 
the mean and standard error of the strain-normalized non-affinity parameter, 4 for PA gel samples 
at various monomer (viz., 7.5% and 15% acrylamide, w/v), and cross-link (between 0.005% and 
0.12% bisacrylamide, w/v) concentrations. The large error bars in the 4 values for the 7.5% acry- 
lamide samples at different bis concentration (standard deviation ~ 38%) arises primarily from 
sample-to- sample variations associated with gels polymerized under (ostensibly) identical exper- 
imental conditions. The error bars for the 15% acrylamide samples are much smaller than the 
7.5% samples, because data in the former case were extracted from a single sample polymer- 
ized at the given acrylamide and bis concentration. Within this relatively large range of values, 
the strain-normalized non-affinity measure does not appear to vary significantly as a function of 
either the density of polymer chains or the network mesh-size, and is evenly distributed around 
the mean 4 calculated over the entire range of PA gels sampled in our experiments. This mean 



value is indicated by the dashed line in Fig. [figure] [4] []4^b). The average values of 4 obtained 
for different monomer concentrations along with their respective standard deviation and standard 



error of each group of measurements are summarized in Table [table] [1][]1 We also used an al- 
ternative approach for calculating the mean ^ in Section 4 of the supporting material available 
online; this alternative approach treated all beads across all samples equally. The results obtained 
by this alternative method were essentially same as the results above. Perhaps not surprisingly, 
we will propose below that this measured non-affinity is largely dominated by inhomogeneities 
formed during synthesis of the PA gels, rather than being dictated by the thermal fluctuations of 
the polymer cross -links. 
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Table 1: Summary of 4 at different acrylamide concentrations. 
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4 Discussion 



4.1 Effectiveness of cross-links 

The measured G' of polyacrylamide gels generally follows predictions of standard theories of 
rubber elasticity, i.e., G' = IvNkgT where /V is the number density of cross-links, kg represents 
the Boltzmann constant, T represents temperature, and v is the efficacy of cross-link. 16 Note the 
additional multiplicative factor of 2, which arises because bisacrylamide is a tetra-functional cross- 
link, v = 1 implies that all cross-links are effective, i.e. the polymer strands attached to each 
cross-link are a part of the homogeneous network. Any unproductive reaction of bisacrylamide 
or inhomogeneity in the network, for example if one of the four polymer strands connected to a 
cross-link is a dangling chain which does not contribute to the elasticity of the polymer network, 
leads to V < iJ 14 l 16 l 22 l Taking the molecular weight of a bisacrylamide as 154, the elasticity of 
a polyacrylamide gel at room temperature can be rewritten as G' = 33.2 x 10 4 vc measured in 



Pascals, where c is the percentage concentration of bisacrylamide. In Fig. [figure] [3] [] 3 b), we see 
that G' is equal to 4.0 x 10 4 c Pa for polyacrylamide gels with 7.5% acrylamide, and 3.0 x 10 5 c Pa 
for 15% acrylamide respectively. Hence v = 0. 12 for the 7.5% and v = 0.9 for the 15% PA gels. 
The higher value of v for 15% acrylamide PA gels is due to the higher polymer chain density in 
this system, suggesting that there is a higher probability for a bis molecule to find an acrylamide 
polymer chain in its neighborhood that would result in an effective cross-link. 

Note though, that this is a simplified scenario; V does not keep increasing indefinitely with 
increasing polymer chain concentration, but levels off as the semi-dilute limit for acrylamide chains 
is reached. The elastic modulus is also strongly affected by the amount of bis present, an excess of 
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which may change the polymer solubility from good to theta solvent and the effective persistence 
length of the acrylamide chains, ■^and may even lead to macroscopic syneresis for sufficiently large 
concentrations of bis. 24 Thus, the relative concentrations of acrylamide and bis may have profound 
effect on the bulk modulus of PA gels, where instead of a linear scaling of the elastic modulus with 



the bis concentration, as seen in our samples (Fig. |[figure] [3] []3f b)), the elastic modulus will level 
off^due to micro-phase separation in the gels. The roles of bis and acrylamide are not limited to 
macroscopic elastic modulus only, but also have significant impact on the microscopic non-affinity 
of PA gels, as we note in the following section. 



4.2 Elastic Inhomogeneities in Polyacrylamide Gels 

The value of 4 characterizes the inhomogeneities in the elasticity of the material. In this section we 
analyze two possible scenarios that give rise to randomness in elastic modulus K (i.e., randomness 
in the shear modulus G' of PA gels), and thus generate non-affine responses in polyacrylamide 
gels. We then compare the predicted non-affinity in these two scenarios with experimental results. 

In the first scenario, the gel is assumed to be nearly ideal. The inhomogeneity is assumed to be 
produced by the intrinsic randomness in the network geometry arising from thermal fluctuations 
frozen into the PA gel at the moment of gelation; in this case, the smallest lengthscale character- 
izing the inhomogeneity would be the mesh size of the network, and the inhomogeneities cannot 
be reduced by improving the synthesis process. In the second scenario, "nonthermal" inhomo- 
geneities are assumed to be introduced during the sample preparation process. For example, this 
effect could arise if the cross-links are not homogeneously distributed; in this second case, elastic 
inhomogeneities could be present at lengthscales much larger than the network mesh size. 

The elastic modulus correlation function in the first scenario can be modeled as 



A K (x) = (8G') z ^8(x), (5) 
where (8G') 2 is the variance of the local shear modulus G', and <^ e is the characteristic mesh size of 
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the network. In this scenario the gelation process is nearly ideal, and the only randomness comes 
from the frozen thermal fluctuations in the liquid at the moment of gelation. Thus the correlation 
of the elastic modulus is characterized by the only length scale in the system, the mesh size £ e , 
which we also take to be the short-distance cutoff of the system because the picture of continuous 
elasticity breaks down below this length scale. Thus the elastic modulus correlations at this scale 
is characterized by a Dirac delta-function in Eq. Q. The Fourier transform of A K is then 

A K (q) = {SG'm (6) 

We plug this correlation function back to Eq. ([6]) to calculate A. To evaluate the integral, one 
has to set a small length scale cutoff, which is as we discussed. Below this length scale the 
polymer network structure cannot be coarse grained, and one cannot characterize the properties 
using continuous elasticity. 

We obtain the non-affine parameter A from the integral (ignoring unimportant 0(1) constant 
pref actor): 

*~(t£)V«?. < 7 > 

The experimentally measured A is plotted as a function of y 2 at two different monomer and 



cross-link concentrations in Fig. [figure] [4] []4^a). Both results are consistent with the theoretical 



prediction that A is proportional to y 2 . However, the magnitude of 4 is of the order of 1 /im 2 . Since 



the mesh size is expected to be of order 10 nm, we obtain 8G' / G' ~ 0(10 ) from Eq. (|7j) and our 
measurements of 4. This value is too large for a nearly ideal "thermal" gel. In a nearly ideal 
gel, the inhomogeneities in network geometry and thus the elasticity come purely from thermal 
fluctuations at the moment of gelation, thus both 8G' and G' are of order kgT times the cross-link 
number density, so one should expect 8G' /G' ~ 0(1). Furthermore, in this scenario A should 
be related to the concentration of cross-links c as A <^ 2 c~ 2 / 3 (because £, e c~ 1//3 ), but this 



behavior is not seen in Fig. [figure] [4] [] 41)), in which 4 is essentially a constant (albeit with a 
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wide scatter). 

Other length scales in the gel will not significantly affect this analysis. The persistence length 
of the polymer chain is even smaller than the mesh size, the small length cutoff of the analysis, and 
thus will not affect the result. The size of the tracer bead, although on the scale of 1 \lm and relevant 



to the problem, only weakly changes the value of A, as we discuss in Section pubsection] [4][4]4.4 
and Section 2 of the supporting material available online. 

The discrepancy between the value of 8G' /G' suggested by the experiment and the theoreti- 
cal value of nearly ideal "thermal" gel suggest that our second scenario may be more realistic for 
these systems. In the second scenario, "nonthermal" inhomogeneities in the distribution of the 
bisacrylamide (cross-link) during the process of polymerization are assumed to exist. These inho- 
mogeneities are frozen in at polymerization and their contribution to the inhomogeneous elasticity 
in the resulting PA gel dominates over the contributions of thermal fluctuations at gelation de- 
scribed in the first scenario, because these "nonthermal" inhomogeneities exhibit greater variance 
and longer correlation length, as we discuss below. These types of heterogeneities have been rec- 
ognized previously in the literature J 23 l 24 l 26 ti HH Briefly, because the hydrophobicity of polymerized 
bisacrylamide is higher than the polyacrylamide chains, the hydrophobic cross-links have a ten- 
dency to aggregate during the sample preparation. ™This effect can generate an inhomogeneous 
spatial distribution of cross-links at length scales longer than the mesh size. 

Thus, in the second scenario, regions with high shear modulus and regions with low shear 
modulus form in the polyacrylamide network as a result of the inhomogeneous distribution of 
cross-links frozen in during the process of sample preparation. We define the length scale to 
characterize the size of this inhomogeneity. The resultant inhomogeneity in the shear modulus may 
then be characterized by a Gaussian correlation function, 

A K (r) = (8G') 2 e (8) 

We then plug this correlation function it back into Eq. Q. The integral is convergent due to the 
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finite range of this correlation function, i.e., so the short length scale cutoff is not needed in this 
scenario. The resulting non-affine parameter is given by 




For a careful derivation of this relation with the exact value of the prefactor, see Section 1 of the 
supporting information available online. 

In fact, considerable effort has been expended over the years to characterize inhomogeneities 
inherent to PA gels. Starting with the pioneering work of Richards and Temple (1971),^ var- 
ious experimental techniques, viz., gel-swelling and permeability studies, ^ 32 ! 33 ! small angle x- 
rayl 23 ! 24 ! 2 ^ and neutron scattering,! 23 * 27 ! quasi-elastic light scattering,! 23 ! 35 ! 38 ! dynamic light scatter- 
ing,! 24 ! 37 ! UV-visible^ and IR spectroscopy,^ NMR spectroscopy, ! 30 ! 31 ! electron micrographs,^ 
have been used to quantify the nature and size of inhomogeneities created in PA gels. Some of 
these ideas have been considered in the context of gel elastic properties,^ as well as under varying 
acrylamide and bis concentration, and different polymerization reaction conditions. The ratio of 
monomer to cross-link concentrations, which determines the relative wettability of acrylamide and 
bis clusters during the polymerization process, as well as the reaction kinetics, all affect the forma- 
tion of dense, heterogeneous clusters of highly cross-linked polymers interspersed with patches of 
sparsely cross-linked polymer chains. The size of these spatial inhomogeneties embedded in the 
more uniform gel matrix has been reported to vary widely from a few nanometers to as much as 
half a micron, with homogeneous regions of comparable length scale in between. 

One may substitute the inhomogeneity correlation length, with the size of the spatial inho- 
mogeneities reported in the aforementioned references. From the literature we find that 5 nm < 
<rG ^ 500 nm, which gives corresponding range of inhomogeneity magnitude of 3 < < 300 
for PA gels over a wide range of monomer and cross-link concentrations. For PA gels synthesized 
under similar preparation conditions as in our experiment, the length scale of inhomogeneities has 
been measured using a nano-indentation method, leading to < 200 nm,^from which we obtain 
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H <7 



4.3 Repeated Shear Measurements 

As part of this study, we explored the effects of cycled measurements on A in the same sample. 
By repeatedly shearing and unshearing a PA gel sample at the same strain, we determined the 
distribution of A for the same set of particles within a single sample. The resultant variation of A 
is not insignificant, though it is considerably less than sample-to- sample error. 

To demonstrate this effect, a PA gel sample is synthesized at 7.5% acrylamide and 0.06% bis 
with 1 \im tracer beads embedded in it. The gel is sheared repeatedly to a strain of 0.2, and A is 
measured each time as shown in Fig. [figure][5][]5| ^a). Error bars reflect the systematic error in our 



measurements. The tracer beads relax back roughly to their original (unsheared) positions once the 
strain is released. The variation in A suggests that some local rearrangement of the polymer net- 
work neighborhood occurs after/during each cycle, perhaps because of the presence of compliant 
chain entanglements or reorganization of the gel-bead interface. These rearrangements permit the 
tracer beads to explore and experience slightly different local environments every time the sam- 
ple undergoes a shear transformation. Non-affinity was slightly different after each shear event. 
The measured standard deviation of 4 (~ 8%) for repetitive shear in the sample is much smaller, 
however, than that measured for different gels prepared under apparently identical experimental 
conditions. 

With respect to non-affinity variation with repeated cyling, we have explored this phenomenon 
under different strains as well as for different polymer gel concentrations. It appears that the 
randomness persists even when a sample gel is sheared repeatedly thirty times. The variation 
in non-affinity parameter appears to be random, independent of the number of times the gel is 
sheared. Chain entanglements, dangling ends, etc. could contribute to this randomness in the 
measured non- affinity,^ and one cannot rule out the possibility that the local environment of the 
tracer micro-beads is subtly distorted due to polymer depletion or adsorption, which might cause 
more/less slippage or sticking of the tracer beads to the surrounding gel matrix under shear.^ED 
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We use this repeated shear technique to calculate the systematic error in our measurements to be 



8% an use this value as the lower bound for all error estimations shown in Fig. |[figure] [4] []4fb). 
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Figure 5: (a) Non-affine parameter, A in a sample PA gel (7.5% acrylamide and 0.06% bis), 
sheared repeatedly under y = 0.2 strain, (b) Average non-affine parameter in a sample PA gel 
measured using fluorescent tracer beads of average diameters of 600 nm, 1 \im and 1.5 [lm. A 
decreases linearly with the size of tracer beads. Measurements shown here were performed on a 
sample PA gel with 7.5% acrylamide and 0.06% bis, sheared ten times at a strain of y = 0.3. (c) 
Elastic shear modulus decreases with increasing initiator and catalyst concentrations, for PA gel 
where the monomer and cross-link concentrations have been kept constant(inset). A decreases 
linearly with increasing initiator and catalyst concentrations. Data are shown here for a 7.5% 
acrylamide and 0.03% bis PA gel. 
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4.4 Tracer Bead-size Dependence 

We also explored the effects of the size of the tracer beads on the magnitude of the non-affine 
parameter, A, using tracer beads of three different sizes, viz., 0.6, 1 and 1.5 /im. The different- 
sized beads are fluorescently labeled such that they are uniquely excited by three different wave- 
lengths of the confocal scanning beam, viz., 488, 568, and 640 nm respectively. We disperse 
these three different-sized beads in a sample PA gel and image them using three different wave- 
length excitation beams in succession during a particular shear event. We see that, for the range 
of bead-sizes used in our experiment, the magnitude of A remains within the range indicated in 



Table [table][l][]l for 7.5% acrylamide PA gels. We also note that though weak, there is a slight 



bead-size dependence of the average value of A measured from repeated shear events. Fig. [fig 



ure][5][]5fo ) plots the average A from eleven repeated shear events at 7= 0.3, for a PA gel at 7.5% 



acrylamide and 0.06% bis. We see that the average A decreases with the diameter of tracer beads. 
A can be fit to a linear function of tracer bead size, with a slope of ~ —0.11 ±0.001/im 2 and an 
intercept of ~ 0.27 ± 0.002. 

Essentially all of the theoretical analysis presented in this paper thus far employed the sim- 
plifying assumption that we can treat the tracer beads as point objects that probe local non-affine 
deformations. However, the size of the tracer bead is comparable to the correlation length of 
the random elastic modulus. In Section 2 of the supporting material available online, we compute 
the corrections due to the finite size R of the bead in a simplified model of electrostatics in random 
media, which a scalar analogy to the elastic problem. In the limit of R — > the non-affine parame- 
ter A smoothly approaches the limit of point probe, while in the limit of R/% ^> I, A approaches 
a different value which is related to the R — ► value by a constant factor of O(l). This simple cal- 



culation is consistent with the experimental observation (Fig. [figure] [5] [] 5 [b)) that A only weakly 



depends on bead-size, R. However the exact dependence is not captured by the calculation. 
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4.5 Effects of Initiator and Catalyst Concentration 

Finally, we explored the effects of reaction kinetics on the strain-normalized non-affinity measure, 
4. To do that, we prepare PA gels with same amount of monomer and cross-link concentration, 
viz., 7.5% acrylamide and 0.03% bis, but with the initiator and catalyst (TEMED and APS, respec- 
tively) concentrations twice and half of the normal amount used. The gel reactions proceed faster 
(slower) as a result, respectively, yielding lower (higher) plateau shear modulus for twice (half) the 
normal initator and catalyst concentrations (inset in Fig. 



[figure] [5] [] 5 c)). The 4 values calcu- 



lated for these samples are still within range of 1.65 ±0.63 jlm 2 , the measured average for 7.5% 
acrylamide PA gels, leading us to believe that measured values of 4 are still dominated by the 
inhomogeneities in these gels. Error bars reflect the standard deviation in the ensempble-averaged 
non-affinity values measured for four scans over each sample volume. Within this prescribed 
range, though, there is a slight inverse dependence of 4 on the concentration of TEMED and APS 



(Fig. [figure] [5] []5 c)) which we do not understand. 



5 Conclusions 

Non-affine deformations under shear are measured in a simple cross-linked gel and are employed 
to provide insight about inhomogeneities in the flexible polymer gels. Results indicate that the non- 
affine parameter, A, which is the mean square non-affine deviation in the PA gels, is proportional to 
the square of the strain applied, in agreement with theoretical predictions based on a linear elastic 
theory treating disorder as perturbations. Interestingly, the magnitude of A is greater than what one 
would expect from theoretical calculations assuming that the PA gels are nearly-ideal and the only 
source of disorder is from the frozen-in thermal fluctuations at gelation. Furthermore, the degree 
of non-affinity appears to be independent of polymer chain density and cross-link concentration. 
Thus, we posit that there are additional built-in inhomogeneities in the PA gels that lead to the large 
non-afflninity we observe. Indeed, there is ample evidence in existing literature of the presence of 
such inhomogeneities in PA gels due to a difference in the hydrophobicities of the bisacrylamide 
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and acrylamide monomers. We also discuss the length scale and magnitude of these additional 
inhomogeneities. Our measurements of non-affinity in PA gels, which are model flexible polymer 
gels, provide a benchmark for the degree of non-affinity in soft materials, and will serve as am 
interesting comparison to non-affinity in more complicated materials such as semi-flexible bio- 
polymer networks. 
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A Non-affine correlation function in inhomogeneous media 

We assume that the gel has a homogeneous bulk modulus, but inhomogeneous shear modulus. 
The inhomogeneity in shear modulus originates in the inhomogeneity of the cross-links due to 
effects such as hydrophobicity, etc. at the moment of gelation, while the monomers are relatively 
homogeneous, and thus the bulk modulus is homogeneous. 

The mean and variance of the shear modulus can be written as 

<G(r)> = G 

(G(r)G(r'))-G 2 = {SG(r)SG(/)) = (5G)VV~, (10) 

where (■ • • ) represent disorder average, £, is the correlation length of the shear modulus, and 8G 
represent the strength of the disorder in G. In this Appendix we drop the ' on shear modulus G' for 
convenience of notation. 

To calculate the non-affine deformation, we use the form of the elastic energy, Eq.(3.25), in 
Ref., 2 which is a first order expansion in small randomness in elastic constant K^i and the resulting 
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small non-affine deformation u 

8U = ^J {Kij U djUidiu k + aji(r)djUidiUi + 2SKij U (r)yjadjUi}, (11) 

where we've only kept first order terms. For isotropic media the elastic constant K^y is related to 
the Lame coefficients as 

Kijki = kdij8 k i + G(8 ik 8ji + 8ii8 jk ) , (12) 

in which G is also the shear modulus, and the bulk modulus is given by B = X + 2G. 
The variation in u to minimize 8V, gives to leading order 

-K ijk idjdiu k (r) = dj8K ijk i(r)y kh (13) 

where the c term vanishes at equilibrium. This equation can be written in momentum space as 

K ijk iqjqiu k {q) = -qj8K ijkl (q)y k i. (14) 

We define the driving force fi(q) = —qj8Kij k i{q)y k i, then the above equation can be written as 

[(X + 2G)V L (q) + GV T (q)] ik q 2 u k (q) =fi(q), (15) 

where V\ k {q) = qiq k /q 2 and V[ k (q) = 8 ik — qiq k /q 2 are the projection operators, and they are or- 
thogonal. Thus the solution of the non-affine deformation as a function of given disorder 8K is 



Ui{q) = 



' -,V L (q) + ^V T (q) 



X+2G ^ G 



\fk(q). (16) 



Using this we can calculate the correlation function of the non-affine deformation Q (r, r') = 
{u m (r)u m (r')) in Fourier space 
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G{q,q') = {u m {q)u m (q')) 



\q jq r {SK ijk i (q) 8k v fxv (q'))7ki7k'i'- (17) 



We assume that there's only disorder in shear modulus G, so the correlation function is (as dis- 
cussed in Appendix A and B in RefJ^J) 

(8Ki jk i(q)8Kj>f k >i>(q')) = (8 ik 8ji + 8u8jk){8i>k>8fi> + 8i>i>8f k >)(8G(q)8G(q')). (18) 

As we already discussed, 



(8G(q)8G(q')) = (27r) 3 5(<7 + </)(5G) 2 (27z;) 3/2 £ V^' 



2ui2 



% \q\ 



(19) 



and suppose the shear deformation has the form (7 is defined as A = / + 7) 



r- 



f 7 xz^ 








(20) 



we arrive at the correlation function 



g(q,q') = (2n) 3 8(q + q')g(q) 



(21) 



with 



G(q) 



4q 2 x f z 
(A + 2G)V 



+ 



<tx + TxTy - 2 TxT z + T y Tz + 4z 



G 2 q 2 



2 k\ 2 



^ z (8G) 2 (2n) 3 / 2 ^e-^. (22) 
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Using this we Fourier transform back to real space and get A as 



A = (8u m (r)du m (r)) 
d 3 q 



(lit) 



rQ{q) 



d 3 q 
(2^ 



, ql + q^-2f x q 2 z +q 2 q 2 z +q 4 z 



(A + 2G)V 



+ 



G 2 q 2 



2|„|2 



e 2 



^ z (8G) 2 (2ny 3/2 % 3 J q 2 smedqddd<j) 



xe 2 



g 2 M 2 j 4(sin0cos0cos(^ 



(A+2G)V 



+ 



GV L 



(sin cos (j)) 4 + (sin cos sin0 sin^)^ 



+(sin sin cos Q) 2 - 2(sin cos cos 0) 2 + (cos (T 



t4(5G) 2 (2^)-^ 2 



1 



167T 1 87T 

+ 



L(A+2G) 2 15 G 2 5 J 



*g(jgg £ 2 | 8 , (gG) 2 2 
5 fc G 2 <> ^i5 fa (A + 2G) 2S 



(23) 



The second term may be small compared to the first term, given that the bulk modulus is large. 

B Corrections for finite-sized bead 

In this experiment, local displacements are probed by beads of size in the order of 1 jlm that are 
embedded in the media, and the displacement of a bead can be different from that of a point probe 
(the position of the point probe is the same as the center of the bead). The correction due to the 
finite size of the bead is of order E,g/R where R is the radius of the bead, proved as follows. 

We simplify the discussion by ignoring the vector nature of the elastic displacements, and 
consider the analogous scalar model that studies a finite sized conductor bead in a inhomogeneous 
electric field. The correspondence of the quantities in electrostatics and elasticity are listed in 



Table [table] [2] [2 147483647] 2 



Let us first consider the simplest case, that a conductor bead of radius R placed at the origin 
r = in a homogeneous external field Ee z , the electric potential 0o outside the bead satisfies the 
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Laplace's equation 



V 2 0o(r)=O, (24) 

and the solution is well known 

0o(r) = Ercos6-E^-J rcos0, (25) 

where 6 is the inclination angle of r measured from e z . The potential of the bead itself is (j)Q = 
fo(r)|r|=* = 0. 

Table 2: Analogy between quantities in electrostatics (scalar problem) and elasticity (vector prob- 
lem) of a finite sized bead placed in an inhomogeneous media in with external field. 



Electrostatics 


Elasticity 


e (permittivity) 
(potential) 
E (field) 


K (elastic modulus) 
u (displacement) 
7 (strain) 



For the case of inhomogeneous media, which has randomness in the permittivity e(r) = Co + 
5e(r) that results in a random potential 0(r) = 0o( r ) + ^0( r )> me Laplace's equation is 



V • [(e + 5e(r))V(^ + 8<j>(r))} = 0. (26) 



This equation can be solved perturbatively, and the zero-th order solution 0o is given in Eq. (25) 
The first order perturbation equation is 



£ V 2 S<Kr) = -(VSe(r)) ■ (V^(r)), (27) 



for the electric potential outside the bead (r = |r| > R), and the boundary condition is that at the 
surface of the conductor bead r = R, the potential 50 (r) = constant. This equation is analogous to 



the first order equation, Eq. ( 15 ) for the elasticity case. This electrostatic equation can be solved 
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using the Green's function defined as 

V 2 Sf(r,r') = --S(r-r'), (28) 
£o 

determining the potential at point r with a point unit charge q = 1 at point r' (outside the bead), 
and the boundary condition is that Sf (r) = constant at r = R. This can be solved using the method 
of image charges, and we obtain the following solution 

v M = J_( + (29) 

4k£q \|r-r'| |r-r"| \r\J 

where q' = R/\r'\ and r" = r'R 2 /\r'\ 2 . We are actually only interested in the potential of the bead 
itself, which is 

#(r') ^(r,r')|| rM 

Ak£q R 

= <30) 

which equals to the potential at the point r = (the center of the bead) if there were no bead placed 
there. 

Then using this Green's function we calculate the first order part of the potential 8§ of the 
bead, 

50 = fi*(r)|| rM = J dr-j^tySeV)) ■ (V0 o (r')), (31) 

Thus the correlation function of the potential 8§ is determined by the correlation function of the 
random permittivity e(r) as 

(S0 2 > = / dv x dr 2 - J^^a M -0o(ri)a 2j -0o(r 2 )^i A;(5e( ri )5e(r 2 )), (32) 

J\ti\>r,\t 2 \>r (4^£oJ z |ri||r 2 | 
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where du = d /dr id . Assuming that (5e(ri)5e(r2)} = (5e) 2 C(ri — r^) where C(v\ — Y2) has the 

Iri-r2l 2 

Gaussian form e 2 ^ 2 as we assumed in the elasticity case. This integral can be evaluated in 
momentum space as 



with 



and 



(^ 2 ) = ^jfj J^3 F ^)F(-^C(d) (33) 



C(q) = J JrC(r) e iqr = (2^ 2 ) 3 /V* 2 l q l 2 / 2 (34) 



F WS [ * +*®f (35) 
J\t\>r 4^|r| 

where Eo = V0 O is the electric field in the zero-th order solution. With some tedious calculation 
we obtain 



F(q) = ^f(qR) (36) 



with 



f(qR) = cosfotf) - 2 ( 6 ~ Q 2R2 ) sin (^) ~ 2< ?*( 6 + Q 2R2 ) cos (^) + ^ V ~ g%g)) (37) 

where g = |q| and Sz is the Sine Integral. Therefore (<50 2 ) can be calculated in spherical coordinate 
as 



where * = <^g. Numerical calculation shows that this integral as a function of R/t, has plateaus of 
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value at small R/t, , at which (Sty 2 ) reduces to the value of point probes, and of value ~ 5.6 

at large R/t,, with slowly varying values in between these two limits. Thus, the correlation of the 
potential is 

{8^)=b { -^E^\ (39) 

where b is is a constant of 0(1) with a weak dependence on R/t,. This is consistent with the 
experimental finding that the non-affine parameter only weakly depends on R/t,, although the 
dependence on R/t, is not correctly captured by the theory. 

C Correlation of non-affine displacements between different 
beads 

The correlation of the non-affine displacements between difference beads decay with distance 
quickly as 1/r, and is thus not visible in the current experiment, as discussed below. 

The real space non-affine correlation function is the Fourier transform of Q{q) as given in 
Eq. (3) in the main paper. Using the correlation function for the inhomogeneous shear modulus as 
given in Eq. (8), again in the main paper, we get 




where Erf is the error function. In the limit r — >• 0, 

G(r)\^o = r 1 ( d §r)\ 2 = A (41) 

reduces to the same point non-affine correlation function A which is the non-affine parameter we 
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use. In the limit r 3> % 




In this experiment, the distance between the beads is greater than the diameter of the bead (1 jlm), 
given that t, is about 0. 1 /lm, we have 

e(r)~y@U<0.L/l, (43) 

and decays as r increases. 

Indeed, for beads separated by 500 nm or more, as in our case, we do not detect any correlations 
in the non-affine displacements between different beads. This result is also in agreement with what 
is observed in the numerical results in RefP 



D Ensemble averaged 4 values for PA gel samples 
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Figure 6: (a) 4 for sample PA gels at 7.5% and 15% acrylamide are plotted at varying bis concen- 
trations. The data points represent the ensemble-averaged values of measurements from several 
different samples prepared, ostensibly, in the same manner. Error bars reflect the fitting error in 



calculating ensemble-averaged g. The dashed line in the main figure indicates the mean of all ^ 



A 



values listed in Table [table] [3] [2 147483647] 3 
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Table 3: List of 4 values for different PA gel samples. Error estimates reflect the uncertainty in 
the linear fits from which 4 are obtained. 



Sample# 


acrylamide conc.[%] 


bis conc.[%] 


fi \m 2 '. 


Error [/lm 2 ] 


1 


7.5 


0.03 


1.18 


±0.08 


2 


7.5 


0.03 


1.92 


±0.96 


3 


7.5 


0.06 


2.08 


±0.15 


4 


7.5 


0.06 


0.61 


±0.28 


5 


7.5 


0.09 


1.15 


±0.58 


6 


7.5 


0.09 


2.03 


±0.94 


7 


7.5 


0.12 


1.69 


±0.50 


8 


7.5 


0.12 


2.72 


±0.25 


9 


7.5 


0.12 


1.50 


±0.41 


10 


15 


0.005 


1.13 


±0.08 


11 


15 


0.01 


2.36 


±0.16 


12 


15 


0.03 


1.60 


±0.23 


13 


15 


0.05 


1.71 


±0.07 



Fig. [figure][6][2147483647]6 a) plots the ensemble average of the strain-normalized non- 



affinity parameter, 4 for different PA gel samples at various monomer (viz., 7.5% and 15% acry- 
lamide, w/v), and cross-link (between 0.005% and 0.12% bisacrylamide, w/v) concentrations. Be- 
cause different PA gel samples have different number of tracer beads dispersed in them, this alter- 
native approach is pursued where every tracer bead across different samples is weighed equally. In 
this method, the mean square non-affine displacement collected from various samples at a given 
acrylamide and bis concentration is plotted against the mean square fitted strain, where all the fitted 
strain values are very close to the externally applied strains. The linear fit of the mean square non- 
affine displacements versus the mean square fitted strains, for different externally applied strains, 



gives the ensemble averaged 4> as shown in Eq. (144b. Error bars in the figure reflect the error in the 



linear fits. Note that the strain-normalized non-affinity parameters obtained from this method are 
very similar to that shown in Fig. 4b in the main paper, and confirm the robustness of our results. 



A 



y ,\r. X* 11 
i=\j=\ 
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i m 

1=1 



(44) 



Here z = 1 , 2, . . . , m labels the PA gel samples at a given acrylamide and bis concentration, with the 
z'-th sample contains Nj beads labeled by j = 1 , 2, . . . ,iVj. 

4 calculated from individual samples at various acrylamide and bis concentrations are listed 



in Table [table] [3][2147483647]3[ along with their respective error estimates. 
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Figure 7: (a) Schematic showing shear deformation of a polymer network. Tracer beads marked 
in gray, green and red indicate the unsheared, affine and non-affine displacement positions under 
shear, respectively. An applied strain, / induces a non-affine displacement, u{r), for a given tracer 
bead in a network, (b) Mean square non-affine displacements, (\u(r)\ 2 ) of tracer beads trapped in 
a polyacrylamide gel scale as the square of the applied strain, y. 
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